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Abstract 

We study the effects of decoherence on some particular bipartite qutrit 
states under the influence of global, collective, local and multilocal de- 
polarizing noise. We show that certain free entangled distillable qutrit 
density matrices become bound entangled or separable and hence convert 
into nondistillable density matrices in global noise. The collective noise 
increases the degree of entanglement of the qutrit bipartite states. Fur- 
thermore, we show that some particular local operation cannot avoid the 
Nondistillability of the distillable states. 
PACS: 03.65.Ud; 03.65.Yz; 03.67.Mn;03.67.Pp 
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1 Introduction 

Entanglement is one of the potential sources of quantum theory. It is the key 
concept and a major resource for quantum communication and computation 
PP. In the last two decades, enormous efforts have been made to investigate 
various aspects of quantum entanglement and its benefits in a number of se- 
tups, such as teleportation of unknown states [2], quantum key distribution 
[2], quantum cryptography [3] and quantum computation [5J |B]. However, a 
realistic quantum system cannot be completely isolated from its environment. 
Generally, the coupling between an environment and a quantum system leads 
to the phenomenon of decoherence and it gives rise to an irreversible transfer 
of information from the system to the environment [Sj [9]. Nevertheless, It 
has also been shown [TU1 HU H2] that the environment can induced entangle- 
ment between two otherwise totally decoupled systems. It is, thus, necessary to 
investigate the behavior of initial state entanglement in the presence of decoher- 
ence. Yu and Eberly [T31 [Tl] showed that entanglement loss occurs in a finite 
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time under the action of pure vacuum noise in a bipartite state of qubits. They 
found that, even though it takes infinite time to complete decoherence locally, 
the entanglement may be lost in finite time. The name "entanglement sudden 
death" (ESD) was tossed for the phenomenon of sudden loss of entanglement. 
The finite time loss of entanglement definitely limits the application of entan- 
gled states in quantum information processing. The phenomenon of ESD is not 
limited only to two qubit entangled state, it is investigated in systems of larger 
spaces such as qutrits and qudits [HJ [TU [T71 HH [HJ [201 HU [22]- A geometric 
interpretation of the effect of ESD is given in Ref. [23] . The experimental ev- 
idences of the phenomenon of ESD have been reported for optical setups [21] 
and atomic ensembles [25] . 

The extraction of a pure entangled singlet state from an ensemble of mixed 
state through local quantum operation and classical communication (LOCC) is 
called distillation [25] . An inseparable state is said to possess some degree of en- 
tanglement. Quantum states are grouped into separable and entangled states for 
qubit-qubit and qubit-qutrit system by using Peres-Horodecki criterion [271 128] . 
According to this criterion, the partial transpose of a separable density matrix 
must has non-negative eigenvalues. The partial transpose of a bipartite den- 
sity matrix p mv<nil over the second qutrit B is given by p^ nv = P mv . nil and 
for the first qutrit, it can similarly be defined. Nevertheless, such a charac- 
terization for higher dimensional bipartite states is difficult [25]. A bipartite 
entangled state can either be free or bound entangled state. A free entangled 
state can be distilled to a singlet state through LOCC whereas a bound en- 
tangled state cannot be, no matter how many copies are available. A bound 
entangled state cannot be used for reliable quantum information processing [30] . 
A nondistillable bound entangled state have positive partial transpose (PPT), 
whereas a negative partial transpose (NPT) states are regarded distillable [50] . 
However, no single criterion is sufficient to detect bound entangled states. For 
example, for a qutrit-qutrit system there are many bound entangled states and 
no single criterion can fully describe all of them [3T]. However, it is shown in 
Ref. [32] that realignment criterion can be used to detect certain bound en- 
tangled states. For a bipartite density matrix p mv , the realignment criterion 
is given by yP) mn v — Pmv n^r A state is separable under realignment crite- 
rion if \\p R \\ < 1 and a PPT state is bound entangled if the quantity \\p R \\ — 1 
is positive. This is important to point out here that realignment criterion is 
not capable to detect all bound entangled states. The phenomenon in which a 
free entangled state becomes nondistillable is called distillability sudden death 
(DSD). It is shown in Refs. [22, 33, 34 that in the presence of dephasing noise 
some free entangled states of qutrit-qutrit system become nondistillable in a fi- 
nite time. The effect of amplitude damping channel on such a system is studied 
in Ref. |35] , where the authors found that a simple local unitary transformation 
can completely avoid DSD. 

In this paper we study the behavior of entanglement by revisiting a par- 
ticular family of density matrices that are considered in Refs. [22j [33] [34] [35] 
of qutrit-qutrit system in the presence of depolarizing noise. We use partial 
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transpose criterion and realignment criterion for our investigation. We consider 
various coupling of the system and environment in which the system is influ- 
enced by global, collective, local or multilocal depolarizing noise. It is shown 
that under certain particular conditions, the depolarizing noise can increase the 
degree of entanglement in the qutrit-qutrit system. For example, under the ac- 
tion of multilocal depolarizing noise, the negativity increases linearly with the 
increase of depolarizing noise of one qutrit's local environment by controlling 
the other qutrit's local environment. On the other hand, the free entangled den- 
sity matrices can become bound entangled or separable and hence nondistillable 
under another particular condition in the depolarizing noise. For example, all 
the free entangled density matrices become bound entangled under the action 
of global environment. Furthermore, it is shown that DSD cannot be avoided 
by performing simple local unitary operation, considered in this paper, for the 
density matrices. Since no definitive criterion for separability or entanglement 
of density matrices of dimensions greater than six is available, we believe that 
our results are valid for the family of density matrices considered in this paper. 



1.1 Qutrit-Qutrit System in a Depolarizing Noise 

We consider a composite system of two qutrits A and B that are coupled to a 
noisy environment both collectively and individually. The qutrits are spatially 
separated and has no direct interaction with each other. The collective coupling 
refers to the situation when both the qutrits are influenced by the same envi- 
ronment and the multilocal coupling describes the situation when each qutrit 
is independently influenced by its own environment. The system is said to be 
coupled to a global environment when it is influenced by both collective and 
multilocal noises at the same time. Let the bases of Hilbert space of each qutrit 
be denoted by |0), |1) and |2). Then the bases of the composite system are given 
in the order |00), |01), |02), |10), |11), |12), |20), |21), |22). 

The dynamics of the composite system in the presence of depolarizing noise 
can best be described in the Kraus operators formalism. The Kraus operators 
for a single qutrit depolarizing noise, which satisfy the completeness relation 
J2i EjEi — I, are given as [35] 




(1) 



with 
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Y=\00l\, Z=\OujO\, (2) 



(3) 



where w = e 2l7r / 3 and p = 1 — e~ rt//2 G [0, 1] is the decoherence parameter. 
The lower and upper limits of p stand, respectively, for undecohered and fully 
decohered cases of the noisy environment. The evolution of the density matrix 
when it is influenced by the global depolarizing noise is given by 

80 8 

P(t) = E E (E? B E?Ef) P (EfEf^Ef^) , (4) 

i=0 j=0 

where — E m eg) 1$, E B =13® E m are the Kraus operators of the multifocal 
coupling of each individual qutrit and Ef B are the Kraus operators of the 
collective coupling that are formed from all the possible combinations of the 
tensor product of the Kraus operators of a single qutrit depolarizing noise in 
the form E m <g) E m . The subscripts m = 0,1, 2,. ..8 stand for a single qutrit 
Kraus operators of a depolarizing noise given in Eq. (fTJ) . The system is said to be 
under the action of global depolarizing noise when both multifocal and collective 
couplings are switched on simultaneously (global = multifocal + collective). We 
consider the following initial density matrix of the bipartite qutrit system 

,„, 2 . . . . a 5 — a , . 

p(0) = -\y + )(y+\ + -<T + + ^-a_ 1 (5) 

where 2 < a < 5. In Eq. ([5]) is a maximally entangled bipartite 

qutrit state given by |*+) = 1/^3 (|00) + 1 1 1 > + |22)), <r + and <j- are sepa- 
rable states, which are given by o+ = 1/3 (|01) (01| + |12)(12| + |20)(20|), ct_ = 
1/3 (| 10) <10| + |21)(21| + |02)(02|). The density matrix of Eq. © is separable 
for 2 < a < 3, bound entangled for 3 < a < 4 and free entangled for 4 < a < 5 
[3"T] . For a further detailed study of the properties of these density matrices we 
refer the readers to Refs. [31] |32j [38] . We will concentrate on the entanglement 
behavior of free entangled density matrix in the range 4 < a < 5. Using initial 
density matrix of Eq. ([5]) in Eq. (0| and taking the partial transpose over the 
second qutrit, it is easy and straightforward to find the eigenvalues. Let the de- 
coherence parameters for multifocal noise of the two qutrits and collective noise 
of the composite system be p\, P2 and p respectively. Then, the eigenvalues of 
the partial transpose of the final density matrix of Eq. ([4]) when only the first 
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qutrit is coupled to the noisy environment are given by 



Ai,3,5 = ( 40- 3pi + \/(41 - 20a + 41a 2 ) (9pi - 8) 2 
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A 2 ,4,6 = ^40-3p 1 -^/(41-20a + 41a 2 )(9p 1 -8) 2 

A 7 , 8l9 = 1 + (6) 

The eigenvalues of the partial transpose of the final density matrix when both 
the qutrits are coupled to multilocal depolarizing noise are given by 

Ai,3,5 = ^(320- 24 P2 + 3 Pl (-8 + 9p 2 ) 

+ V(41 - 20a + 4a 2 )(-8 + 9pl ) 2 (-8 + 9p2) 2 ) 

A 2 ,4, 6 = ^(320- 24 P2 + 3 Pl (-8 + 9p 2 ) 

- V(41 - 20a + 4a 2 )(-8 + 9pi) 2 (-8 + 9p 2 ) 2 ) 

A 7 ,s,9 = -^{l2S + 2Ap 1 + 2Ap 2 -27p 1 p 2 ) (7) 

It is easy to see that Eq. ([7]) reduces to Eq. ^ , if we switch off the coupling 
of second qutrit with its local environment. The eigenvalues when both qutrits 
are coupled to the same environment, that is, for collective coupling of the 
environment and system are obtained by replacing pi = p 2 — p in Eq. ([7]). The 
only eigenvalues that possibly become negative in Eqs. (j6|), © and hence for 
the collective noise are those which are given by A 2j 4 ; 6- To observe the behavior 
of entanglement and sec whether the phenomena of DSD and ESD happen in 
these cases, we use the partial transpose criterion and realignment criterion in 
the following. 

The degree of entanglement in NPT states is quantified by using negativity 
[27} [28] . It is given by the sum of the absolute value of the negative eigenvalues 
of the partial transpose of a density matrix. The negativity for local, multilocal 
and collective noise becomes 

N{ P ) = max{0, |A 2 , 4 , 6 |} = max{0, |3A 2 |} . (8) 

The negativity for different values of the parameter a, when only one qutrit is 
coupled to the environment, is plotted against p\ in Fig. 1. It can be seen that 
the negativity becomes zero at a finite value of p\ (around p\ = 0.2) only for 
density matrices that correspond to the upper limit of the parameter a and hence 
becomes PPT that might possess some degree of entanglement. For density 
matrices in the lower limit of a, the negativity first decreases to a minimum 
value, however, still positive and then increases with pi, thereby increasing 
the degree of NPT (free entanglement). The partial transpose criterion fails 
to detect the behavior of entanglement in density matrices of large a beyond 
certain values of p\ . Further investigation of the entanglement in these states 
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Figure 1: (color online) The negativity for various values of parameter a is 
plotted against the decoherence parameter p\ , for the case when only one qutrit 
undergoes decoherence. The realignment criterion \\p || — 1 for single qutrit 
decoherence is also plotted for a = 5. 

can be carried out by using realignment criterion. The quantity \\p R \\ — 1 for 
density matrices of large a is also plotted in Fig. 1. The positive value of this 
quantity in the range 0.21 < pi < 0.288 shows that these density matrices are 
bound entangled. However for p\ > 0.288, the realignment criterion also fails 
to detect the possible entanglement for these states. 

The negativity for the case when both qutrits are individually coupled to its 
own environment, that is, for multilocal coupling is plotted in Fig. 2 against 
the decoherence parameter pi for decoherence parameter p\ = 0.1. It can be 
seen that neither DSD nor ESD occurs in any density matrix (4 < a < 5) as 
the negativity is positive for all of them. For density matrices of lower values 
of a in the range 4 < a < 5, the negativity increases linearly with p2 and 
hence the degree of entanglement. On the other hand, for density matrices that 
correspond to intermediate and upper values of a, it first decreases, however 
remain positive, and then increases linearly. It can be shown that the negativity 
increases linearly for all density matrices, irrespective of the value of a (4 < 
a < 5), for larger values of the decoherence parameter p\. Thus, both DSD 
and ESD can be completely avoided in the case of multilocal coupling of the 
system and environment by controlling the environment of one qutrit at least. 
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Figure 2: (color online) The negativity for various values of parameter a is 
plotted against the decoherence parameter p2 when both the qutrits are coupled 
to their local environments. The value of decoherence parameter p\ = 0.1. 



Nevertheless, for pi < 0.1, some density matrices that correspond to large values 
of a become PPT. For example, for p\ — 0.05, the negativity of density matrix 
that corresponds to a — 5 becomes zero at pi — 0.165. The negativity and the 
realignment criterion for such density matrix is plotted in Fig. 3. It shows that 
the density matrix becomes bound entangled for 0.165 <pi < 1. 

In Fig. 4, we have plotted the negativity for various density matrices against 
the decoherence parameter p for the case of collective noise. As can be seen from 
the figure, non of the density matrices undergoes DSD or ESD. According to 
partial transpose criterion, all the density matrices in the range 4 < a < 5 
are free entangled under the action of collective depolarizing noise. Though the 
negativity for each density matrix drops to a minimum, but positive value, in the 
range of lower values of decoherence parameter. The degree of NPT increases 
for all density matrices with the values of decoherence parameter p in the large 
values limit. 

The general form of eigenvalues for the global noise have quite lengthy ex- 
pressions, instead of writing their expression, we prefer to see their behavior by 
plotting them against decoherence parameter p. Such a plot for p\ = p2 =0.5 
and a — 4.3 is shown in Fig. 5. Since all the eigenvalues for the whole range of 
decoherence parameter are positive, the density matrices under global noise are 
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Figure 3: (color online) The negativity and realignment criterion \\p R \ \ — 1 for 
the multifocal coupling are plotted against the decoherence parameter p2- The 
values of other parameters are set to p\ = 0.05, a = 5. 

PPT. The realignment criterion also fails to detect any possible entanglement 
as the quantity 1 1 1 1 — 1 for the range of interest of the parameter a against the 
decoherence parameter p is negative. We conclude that the distillable free en- 
tangled qutrit-qutrit states becomes completely nondistillable in the presence of 
global depolarizing noise. Whereas under the action of global dephasing noise, 
it is shown 34j that these density matrices are free entangled in certain limited 
time and undergo DSD at a later time. 

To see the effect of local unitary operation on the dynamics of the state 
of Eq. O we consider the unitary operator U — I3 ® 9, with 9 = |0)(1| + 
|1)(0| + |2)(2|. When this operator is applied locally to the state of Eq. O 
it converts the maximally entangled state |^+) into another maximally en- 
tangled state given by |*+) = 1/73(100) + |11) + |22)) and the two separa- 
ble states (j _|_ , er_ into other two separable states, respectively, given by rr + = 
1/3 (|00)(00| + |12)(12| + |21)(21|), = 1/3 (|11)(11| + |20)(20| + |02)(02|). The 
density matrix after U acts on p(0) becomes 



The evolution of ag in the noisy environment and its partial transpose can 



a e - Up(0)Ui = |*+)(*+| + ^+ + 




(9) 
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Figure 4: (color online) The negativity for various values of a is plotted against 
decoherence parameter p when the system evolves under collective depolarizing 
noise. 

be found straightforwardly in the same manner as done for the state p(t). It 
is found that all the eigenvalues for each individual case, considered for the 
coupling of p(t), remains unchanged. In conclusion, the local operation that we 
consider here does not change the behavior of entanglement in the presence of 
any of the aforementioned coupling under the action of depolarizing noise. It 
is important to state that under the action of amplitude damping noise, a local 
unitary operation changes the behavior of entanglement and completely avoids 
DSD [35]. 

2 Summary 

We study the dynamics of entanglement for particular bipartite qutrit density 
matrices under global, collective, multilocal and local depolarizing noise. We 
show that unlike the cases of dephasing |34j and amplitude damping noises |35j , 
the influence of depolarizing noise is completely different. Using partial trans- 
pose criterion and realignment criterion, it is shown that only those density 
matrices that correspond to large values of a becomes PPT and thus nondistill- 
able in the local depolarizing noise. In the case of multilocal depolarizing noise, 
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Figure 5: (color online) The eigenvalues for the case of global noise are plotted 
against the decoherence parameter p. The values of other parameters are set to 
Pi = P2 = 0.5, a = 4.3. 

both ESD and DSD can be avoided by controlling one or the other local envi- 
ronment. We also show that instead of ESD and DSD to occur, the degree of 
entanglement for free entangled density matrices under collective depolarizing 
noise increases with the increasing value of decoherence parameter. All the free 
entangled density matrices in the specified range of parameter a becomes PPT 
under the action of global depolarizing noise. In conclusion, we show that the 
free entangled distillable density matrices convert into bound entangled or sep- 
arable density matrices and thus become completely nondistillable under global 
noise. Also, the depolarizing noise can be used to increase the degree of entan- 
glement in the specific family of density matrices that correspond to a particular 
values of the parameter a under certain couplings with the environment. Fur- 
thermore, It is shown that the local operation, we used in this paper, does not 
change the dynamics of entanglement under each coupling of the system and 
environment for the case of the specific family of the density matrices. Since 
there is no single known definitive criterion for separability and entanglement 
in states of dimension greater than six, therefore, we emphasis that the validity 
of our results is true for the family of density matrices that we consider here 
and may not be generalized to all bipartite qutrit states. 
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Figure Captions 

Figure 1. The negativity for various values of parameter a is plotted against 
the decoherence parameter p\ , for the case when only one qutrit undergoes de- 
coherence. The realignment criterion \\p R \ \ — 1 for single qutrit decoherence is 
also plotted for a — 5. 

Figure 2. The negativity for various values of parameter a is plotted against 
the decoherence parameter p 2 when both the qutrits are coupled to their local 
environments. In the inset of the figure the realignment criterion is plotted 
against decoherence parameter p 2 , which shows that the density matrices be- 
come bound entangled for p^ = 0.165 < p 2 < 1. 

Figure 3. The negativity and realignment criterion \ \p R \ \ — 1 for the multilocal 
coupling are plotted against the decoherence parameter p 2 ■ The values of other 
parameters are set to p\ = 0.05, a — 5. 

Figure 4. The negativity for various values of a is plotted against decoherence 
parameter p when the system evolves under collective depolarizing noise. 
Figure 5. The eigenvalues for the case global noise are plotted against the deco- 
herence parameter p. The values of other parameters are set to p\ = p 2 = 0.5, 
a = 4.3. 
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